Evolution of Surface Deformations of Weakly-Bound Nuclei in the Continuum 
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We study weakly-bound deformed nuclei based on coordinate-space Skyrme Hartree-Fock- 
Bogoliubov approach , in which a large box is employed for treating the continuum and surface 
diffuseness. Approaching the limit of core-halo deformation decoupling, calculations found an ex- 
otic "egg"-like structure consisting of a spherical core plus a prolate halo in '^*Ne, in which the 
resonant continuum plays an essential role. Generally the halo probability and the decoupling effect 
in heavy nuclei are reduced compared to light nuclei, due to denser level densities around Fermi 
surfaces. However, deformed halos in medium-mass nuclei are possible with sparse levels of negative 
parity, for example, in ^^'^Ge. The surface deformations of pairing density distributions are also 
influenced by the decoupling effect and are sensitive to the effective pairing Hamiltonian. 

PACS numbers: 21.10.Gv, 21.10.Pc, 21.60.Jz 



In exploring the limits of nuclear landscape, exotic phe- 
nomena such as halo states and associated soft mode 
excitations are expected near the particle drip- lines 
The halo structure occurs as quantum many-body sys- 
tems approaching decay thresholds. It is characterized 
by tremendous surface diffuseness and has been observed 
in weakly-bound nuclei and molecules [3]. In particular, 
drip-line nuclei are weakly-bound superfluid systems, in 
which the pairing induced continuum coupling plays an 
essential role [1, |3| • Due to the dilute density in the halo 
region, the BEC-like di-neuron pairing is expected to be 
coexisted with the BCS pairing ^5', 6]. A crucial issue for 
theoretical descriptions of drip-line nuclei is to be able 
to self-consistently treat continuum effects and density 
diffuseness. 

There have been many studies of drip-line nuclei with 
assumed spherical shapes and have brought plentiful new 
insights, for example, the core- halo decoupling pair- 
ing anti-halo effect \A, shell quenching |8l] , and contin- 
uum coupling effect [1, 0, H, [l^l ■ Beyond the spherical 
shape, deformed nuclei can provide unique chances to 
study many-body correlations through breaking geomet- 
ric symmetries, e.g., shape evolutions and associated col- 
lective motions. Especially in drip-line nuclei, intriguing 
surface deformations of halo structures could happen due 
to the core-halo deformation decoupling ll|, |l2| . Another 
interesting phenomena is the expected deformation dif- 
ference between neutrons and protons near the neutron 
drip line, which implies unusual isovcctor quadrupole col- 
lective modes [ll|. Besides, in the weak-coupling region, 
the deformation of pairing densities (or anomalous den- 
sity) could be different from particle densities, while such 
properties can be studied by pair transfer experiments 
and have been rarely discussed. 

To describe weakly-bound nuclei, especially for 
medium-mass and heavier nuclei, Hartree-Fock- 
Bogoliubov (HFB) theory is the theoretical tool of 
choice by properly taking into the continuum effect. 



For deformed nuclei in the continuum, however, exact 
HFB solutions with outgoing boundary conditions are 
too difficult. The Green function method has made 
remarkable progress in solving the deformed HFB 
equation 13|, while self-consistent calculations are still 
missing. On the other hand, the coordinate-space HFB 
approach, referred to as the discretization method, 
has been demonstrated to be able to treat the resonant 
and non-resonant continuum rather accurately 
compared to the exact Gamow HFB solutions 
For weakly-bound deformed nuclei 
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calculations have 
to employ a large coordinate-space box that is crucial 
for describing the spatial extension and the continuum 
discretization 0, Recently such very expensive 

calculations of axial symmetric nuclei have been re- 
alized by performing hybrid parallel calculations on 
supercomputing facilities 16|. Since the computing 



capabilities and resources are developing so fast that the 
coordinate-space HFB approach is taking an advantage 
to study deformed weakly-bound nuclei. 

According to systematic calculations of the nuclear 
landscape, there are several deformed regions close to 
drip lines [171 118| . Experiments have been able to reach 
the neutron drip line in the light mass region to Mg 
and Al isotopes [l9j |. which are predicted to be de- 
formed 
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Extensive studies have been done to 
look for deformed halo structures in light nuclei experi- 
mentally [13, [in and theoretically [H [H [H]. For the 
halo mechanism in weakly-bound heavy nuclei, theoret- 
ical studies are mostly done with spherical shapes and 
are model dependent |24l428l |. Actually, the on-going ex- 
perimental facilities such as FRIB are expected to reach 
80% of all possible isotopes for Z <90 [29|, which will 
provide an unprecedent opportunity for exploring heavy 
drip-line nuclei. In this context, we systematically stud- 
ied the even-even weakly-bound deformed nuclei from 
light to heavy mass region based on the self-consistent 
coordinate-space Skyrme HFB approach. While the odd- 



N halo nuclei with much reduced pairing correlations and 
continuum couplings due to one quasiparticle excitation 
are not involved here. 

In this work, the Sk yrm e-HFB equation is solved by the 
HFB-AX code 3 [3, ll^ within a 2D lattice box, based 
on B- spline techniques for axial symmetric deformed nu- 
clei 3l|. To obtain sufficient accuracy, the adopted 2D 
box size is 30x30 fm. The maxima mesh size is 0.6 fm 
and the order of B-splines is 12. This is the first deformed 
HFB calculations with such a large box size, while in lit- 
eratures the adopted 2D box sizes are usually less than 
20 fm. From 20 fm to 30 fm, the estimated computing 
cost will be increased by 40 times [32] • For calculations 



employing large boxes and small lattice spacings, the dis- 
cretized continuum spectra would be very dense and pro- 
vide good resolutions. Because the computing cost is 
extremely high, the hybrid MPI-|-OpenMP parallel pro- 
gramming is implemented to get converged results within 
a reasonable time 
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Calculations were performed in 
China's top supercomputer Tianhe-IA. For the particle- 
hole interaction channel, the SLy4 force [s^ is adopted as 
it is one of the mostly used parameterizations for neutron 
rich nuclei. For the particle-particle channel, the density 
dependent pairing interactions with the mixed and sur- 
face variants are used 3^. The pairing strengthes are 
fitted to the neutron gap of ^^°Sn. 




FIG. 1: (Color online) The density profiles of protons and neu- 
trons in ^**Ne, ""'''''Mg, ^^°Ge, and ^^"Er, which are obtained 
from calculations with the surface pairing. The density dis- 
tributions are shown along the cylindrical coordinates z-axis 
and r-axis, respectively. The dashed lines show the densities 
contributed by near-threshold qp resonances as denoted in 
FigH 



We have searched for deformed halo structures along 
the neutron drip line around three deformed regions: (I) 
the Ne and Mg isotopes, (II) the medium mass region 
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FIG. 2: (Color online) The two-dimensional distributions of 
the neutron density (a) and neutron pairing density (b) in 
^*Ne, which are displayed on the logarithm scale. The shapes 
of the core and halo are denoted by dashed lines. 



around Z=30, and (III) the heavy mass region around 
Z—60. In Figlll to show the surface deformations, the 
density profiles of ^^Ne, ^^Mg, ""Ge, and ^^^Ei have 
been displayed along the cylindrical coordinates z-axis 
(the axial symmetric axis) and r-axis, respectively. These 
four drip-line nuclei all have considerable pairing fields. 
Actually ^^°Ge and ^^^Er are calculated to have the most 
significant halo characteristics around Z=,iO and Z=60, 
respectively. We have to note that predictions of the neu- 
tron drip line can be dependent on the effective Skyrme 
interactions and pairing interactions fl7] . Results shown 
in Fig. [1] are obtained with the surface pairing since halo 
phenomena are more pronounced than with the mixed 
pairing. In Fig. [U '^^Ne has the most pronounced halo 
structure. Generally halo structures gradually fade away 
as nuclei towards heavier mass region. It can be seen 
that Mg has a prolate core and a prolate halo, while 
in relativistic Hartree-Bogoliubov calculations it has a 
prolate core and a slightly oblate halo 12|. In ^^Ne 
is very surprising to see an exotic "egg' 
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-like structure, 
in which the neutron halo distribution is well prolately 
deformed but the core is spherical, as shown in FiglJJa). 
Inside the core with a radius of less than 10 fm, both the 
proton and neutron density distributions have a spheri- 
cal shape. Such an exotic structure has never been ob- 
tained in earlier self-consistent calculations, emphasizing 
the importance of precisely treating the subtle interplay 
between continuum couplings and surface deformations 
in our approach. The "egg" -like structure is also a mani- 
festation of the completely core-halo deformation decou- 
pling effect. In Figl2jb), the neutron pairing density dis- 
tribution of -^^Ne also shows the "egg" -like structure, re- 
flecting the decoupling effect. As nuclei become heavier, 
however, one can see that surface deformations are more 
similar to core deformations, indicating the weakening of 
the core-halo deformation decoupling effect. The most 
significant halo structure in "^^Ne is also consistent with 
an earlier point that halos can be hindered by deformed 
cores [sHi- To study the halo mechanism, in Fig. [1] we also 
displayed the densities contributed by the quasiparticle 
(qp) resonances close to Fermi surfaces. It can be seen 
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FIG. 3: (Color online) The smoothed occupation numbers of r2'^=l/2^ neutron qp states of ^*Ne, '^''Mg, Ge and Er isotopes. 
The shadow areas indicate the near-threshold qp resonances. The neutron Fermi surface A„ are also given. The arrows in Ge 
and Er panels are given for guiding eyes. 



that qp resonances near thresholds are mainly responsible 
for the halo occurrences and surface deformations. In our 
calculations, the halo phenomena are more pronounced 
with the surface pairing than with the mixed pairing. 
This can be understood that surface pairing correlations 
and surface diffuseness can be enhanced due to contin- 
uum couplings With the mixed pairing, however, the 
continuum coupling effect can not be competitive to the 
pairing anti-halo effect. 

To understand the deformed halo mechanism in per- 
spective of qp spectra, we have displayed the qp spectrum 
with smoothed occupation numbers in Fig|3l It is known 
that the discretized qp resonance has roughly the Breit- 
Wigner shape For demonstration, we smoothed the 
qp occupations vf with a Lorentz shape function and 
a smoothing parameter of 0.2 MeV (the smoothing is 
described in Ref. In Fig-3, the smoothed occupa- 

tions of fl'"=l/2^ neutron qp states of '^^Ne and ^''Mg are 
shown, while contributions from high f2 states are much 
suppressed. It can been seen that the main difference 
between "^^Ne and ''^Mg is that '^^Ne has a stronger qp 
resonance near the Fermi surface. Although it is difficult 
to discuss the details of such qp resonances due to angular 
momentum mixing, qualitatively one can infer that the 
deformed halo in ^^Ne can be ascribed to a p-wave con- 
figuration, based on the nodes in density distributions 
as shown in FigH] It has been pointed out that nega- 
tive parity states are responsible for deformed halos [llj . 
However, in ''^Mg a 1/2^ resonance from 2^3/2 around 
3 MeV is too high to contribute to the halo. Compared 
to '^^Ne, as shown in Figdl ''"Mg has the same magic 
neutron number of N=28, but halo disappears because 
pairing correlations are almost gone in our calculations 



and thus near-threshold resonant continuum vanishes. It 
can be concluded that the spatial extension of deformed 
halo structures are closely related to the magnitudes of 
near-threshold qp resonances. 

In Figl21 the systematics of the smoothed qp spectrum 
of Ge and Er isotopes are also displayed to understand 
the development of deformed halos in heavy nuclei. One 
can see that as isotopes close the neutron drip line, the 
weakly-bound qp states move away from the Fermi sur- 
face collectively. Simultaneously, small qp resonances 
gradually develop and decouple from cores (or bound 
states), resulting in deformed halo structures. It can 
be understood that a large phase space between weakly- 
bound states and the Fermi surface is essential for the de- 
velopment of near-threshold small qp resonances. Com- 
pared to ^^Ne and ^'*Mg, level densities around Fermi 
surfaces become larger in heavy nuclei. Therefore, such 
phase spaces can not be very large and halos would be 
less significant. This is consistent with the fact that 
heavy nuclei have larger surface tension due to denser 
level densities around Fermi surfaces and thus have less 
halo probabilities. However, ^^°Ge actually has a simi- 
lar halo asymptotic to '''^Mg although it is much heav- 
ier. This is because the level density of negative par- 
ity 1/2^ states is smaller in ^^°Ge compared to ^''Mg, 
as shown in Fig. [31 Indeed, the l/2~ state around 3.5 
MeV from lft.11/2 (the only negative subshell from A^=50 
to A^=82) which doesn't contribute to deformed halos. 
In contrast, the deformed halo is mainly contributed by 
near-threshold small l/2~ qp resonances with compo- 
nents from the intruder 2/7/2. Due to the phase space 
decoupling as shown in Fig. [3] and continuum couplings, 
the magnitudes of near-threshold qp resonances and then 
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FIG. 4: (Color online) Quardpole deformations of neutron 
densities (dashed line), neutron pairing densities (solid line), 
proton pairing densities (short-dashed line), and deformation 
differences (isovector deformation) between neutrons and pro- 
tons (thick solid line). Calculations of Ge and Er isotopes are 
done with the surface pairing and mixed pairing interactions 
respectively. 



the spatial extension are enlarged. This halo mechanism 
is beyond the conventional s or p-wave halo. Based on 
the FRIB's projection jl^, the neutron drip line of Ge 
could be reached experimentally and hopefully deformed 
halos around this region can be discovered. 

Towards the neutron drip line, the deformation differ- 
ence between neutron and proton density, i.e., the isovec- 
tor deformation, is an interesting topic and has been 
studied systematically by the HFB approach based on 
HO basis It is known that generally neutron distri- 
butions are less deformed compared to protons towards 
the neutron drip line, i.e., the isovector deformations are 
negative, implying interesting isovector quadrupole col- 
lective modes [11]. In this work, we have calculated 
the quadrupole deformations /32 of both particle den- 
sity and pairing density distributions using the formula 
>/< > in Ref. 



^2 = (47r/5)< r^Fao >/< > in Ref. [171. The isovec- 
tor deformations, deformations of neutrons and pairing 
densities are displayed in FigH] One can see that similar 
negative isovector deformations are obtained in calcula- 
tions with both the mixed and surface pairing interac- 
tions. In particular, the absolute isovector deformations 
tend to increase towards the neutron drip line with the 
mixed pairing. The neutron drip-line of Ge and Er ob- 
tained by the surface pairing has two more neutrons in 
addition to the mixed pairing. One can see that defor- 
mations of neutron pairing density distributions obtained 
in surface pairing calculations increase near the neutron 



drip line, in contrary to the case of mixed pairing. The 
deformations of neutrons are also slightly influenced with 
the surface pairing. With the surface pairing for '^^Ne and 
''^Mg, the isovector quadrupole deformations are 0.24 and 
—0.09, respectively, which are larger than that of ^^''Ge 
and ^^''Er. The corresponding deformations of neutrons 
and neutron pair densities are 0.24(0.48) and 0.25(0.36), 
respectively. This is also shown in FigH] that neutron 
pairing density has a larger deformation and spatial ex- 
tension than neutron density, for which the non-resonant 
continuum may play a role [3, [2^ . It has to be noted that 
in the "egg" -like structure the isovector deformation has 
a unusually large positive value. In general the deforma- 
tions of neutron pairing densities are different from that 
of neutron densities and are too sensitive to the choice of 
effective pairing interactions to give specific conclusions. 

In summary, we have studied the weakly-bound de- 
formed nuclei based on the self-consistent coordinate- 
space Skyrme Hartree-Fock Bogoliubov approach, to look 
for new insights by precisely solving the continuum with 
discretization within a large box about 30 fm. We 
found an exotic "egg" -like structure for the first time 
with a spherical core plus a prolate halo, for which the 
near-threshold resonant continuum contributions are es- 
sential. The halo probability and core-halo deformation 
decoupling effect have been studied systematically from 
light to heavy nuclei, which turned out to be depend very 
much on the level density oift^ = 1 /2~ states near Fermi 
surfaces. The obtained large isovector deformations im- 
ply unusual isovector quadrupole excitation modes and 
further experiment observations will be very important. 
The core-halo deformation decoupling effect is also pre- 
sented in the pairing density distributions. However, the 
deformations of pairing density distributions are different 
from particle densities and are sensitive to the effective 
pairing Hamiltonian. 
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